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Abstract 

We consider the mixed boundary value problem, or Zaremba's problem for 
the Laplacian in a bounded Lipschitz domain in R", n > 2. We decompose 
the boundary = D U N with D and N disjoint. The boundary between D 
and is assumed to be a Lipschitz surface in dQ. We specify Dirichlet data 
on D in the Sobolev space W^''^{D) and Neumann data in U\N). Under these 
conditions, we find go > 1 so that the mixed problem has a unique solution 
and the non-tangential maximal function of the gradient lies in U'[dVt) for 
1 < p < qq. We also obtain results when the data comes from Hardy spaces for 
p = 1. 
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1 Introduction 

Over the past thirty years, there has been a great deal of interest in studying 
boundary value problems for the Laplacian in Lipschitz domains. A fundamen- 
tal paper of Dahlberg [8] treated the Dirichlet problem. Jerison and Kenig [15] 
treated the Neumann problem and provided a regularity result for the Dirichlet 
problem. Another boundary value problem of interest is the mixed problem or 
Zaremba's problem where we specify Dirichlet data on part of the boundary 
and Neumann data on the remainder of the boundary. To state this boundary 
value problem, we let Q he a bounded open set in R" and suppose that we 
have written d^l = D D N where D is an open subset of the boundary and 
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N = dU \ D. We consider the boundary value problem 



r An = 0, 



in n 
on D 



< 



u = fn, 
du 



on iV. 



(1.1) 




The study of the mixed problem in Lipschitz domains appears as problem 3.2.15 
in Kenig's CBMS lecture notes |16j . Recall that simple examples show that we 
cannot expect to find solutions with the gradient in of the boundary. For 
example, the function Re y/z on the upper half-plane has zero Neumann data 
on the positive real axis and zero Dirichlet data on the negative real axis but 
the gradient is not locally square integrable on the boundary of the upper half- 
space. This appears to present a technical problem as the standard technique for 
studying boundary value problems has been the Rellich identity which produces 
estimates in L^. 

In 1994, one of the authors observed that the Rellich identity could be 
used to study the mixed problem in a restricted class of Lipschitz domains. 
Based on this work and the methods used by Dahlberg and Kenig to study 
the Neumann problem [9], J. Sykes |30l l3T] established results for the mixed 
problem in Lipschitz graph domains. L Mitrea and M. Mitrea [23] have studied 
the mixed problem for the Laplacian with data taken from a large family of 
function spaces. Results have also been obtained for the Lame system [3] and 
the Stokes system [4J. More recently, Lanzani, Capogna and Brown [TH] used a 
variant of the Rellich identity to establish an estimate for the mixed problem in 
two dimensional graph domains when the data comes from weighted spaces 
and the Lipschitz constant is less than one. The present work also relies on 
weighted estimates, but uses a simpler, more flexible approach that applies to 
all Lipschitz domains. 

Several other authors have treated the mixed problem in various settings. 
Verchota and Venouziou [32] treat a large class of three dimensional polyhedral 
domains under the condition that the Neumann and Dirichlet faces meet at an 
angle of less than vr. Maz'ya and Rossman [20l Ell [22j have studied the Stokes 
system in polyhedral domains. Finally, we note that Savare [28] has shown 
that solutions may be found in the Besov space -Bg^*^ in smooth domains. This 
result seems to be very close to optimal. The example Re ^/z described above 
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shows that we cannot hope to obtain an estimate in the Besov space -63^2- 

We outline the rest of the paper and describe the main tools of the proof. 
Our first main result is an existence result for the mixed problem when the 
Neumann data is an atom for a Hardy space. We begin with the weak solution 
and use Jerison and Kenig's results for the Dirichlet problem and Neumann 
problem [15j to obtain estimates for the gradient of the solution on the interior 
of D or N. This leads to a weighted estimate where the weight is a power of the 
distance to the common boundary between D and A^. The estimate involves a 
term in the interior of the domain. We handle this term by showing that the 
gradient of a weak solution lies in U'[^) for some p > 2. The LP{Q) estimates 
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for the gradient of a weak solution are proven in section [3] using the reverse 
Holder technique of Gehring and Giaquinta. Using this weighted estimate for 
solutions of the mixed problem, we obtain existence for solutions with Hardy 
space data by extending the methods of Dahlberg and Kenig Uniqueness 
of solutions is proven in section [5l 

With the Hardy space results in hand, we establish the existence of solutions 
to the mixed problem when the Neumann data is in L^{N) and the Dirichlet 
data is in the Sobolev space W'^''p{D). This is done in sections [6] and [7] by 
adapting the reverse Holder technique used by Shen to study boundary value 
problems for elliptic systems [29]. The novel feature in our work is that we are 
able to substitute an estimate in Hardy spaces while Shen's work begins with 
existence in L?. 

2 Definitions and preliminaries 

We will consider several formulations of the mixed problem (jl.ip . Our goal 
is to obtain solutions where the gradient lies U'{dVt) for 1 < p < qq for some 
go > 1- However, our argument begins with a weak solution where the gradient 
lies in L^(r2). We will show that under appropriate assumptions on the data, 
this solution will have a gradient in U'{d^). 

We describe a weak formulation of the boundary value problem (jl.ip . The 
results of section [3] will hold for solutions of divergence form operators. Thus, 
we define weak solutions in this more general setting. For k = 1,2, ... , we 
use W^^P{Q) to denote the Sobolev space of functions having k derivatives in 
LP{Vl). For D a subset of the boundary, we let w]f'{Vl) be the closure in 
M^^'^(r2) of smooth functions which are zero on D. In the weak formulation, we 
take the Dirichlet data fi) from the Sobolev space W^^'^{Vt). We let W]{^''^{d9) 
be the restrictions to dVl of the space W^' (O). We define ' {dVl) to be 
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the dual of Wj^ ' (917). The Neumann data /at will be taken from the space 
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W~^'^'\d^). liA{x) is a symmetric matrix with bounded, measurable entries 
and satisfies the ellipticity condition M|^p > A{x)^ ■ £^ > M~^|^p for some 
M > and all £ R", we consider the problem 

divAVu = 0, in O 

u = /d, on D 

AVu ■ V = /at, on N. 

1 2 

We say that u is a weak solution of this problem if u — fi) G Wjj (Q) and we 
have 

/ AVu-Vvda = -{fN,v)9n, veW'/in). 
Jn 

To establish existence of weak solutions of the mixed problem, we will assume 
the coercivity condition 

MIlHu) < c\\Vu\\L2^n), n G (2.1) 
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Under this assumption, the existence and uniqueness of weak solutions to (jl.ip 
is a consequence of the Lax-Milgram theorem. In our appUcations, 17 wih be a 
connected, bounded domain whose boundary is locahy the graph of a Lipschitz 
function and D will be an open subset of the boundary. These assumptions are 
sufficient to ensure that (12. ip holds. 

If /at is a function on A^, then we may identify /tv with an element of the 
space W^^^'^''^{d^) by 

{fNA)du = [ fN^da, 4> G w]i^'\dn). 

JN 

1 /2 2 

From Sobolev embedding we have Wjj ' {do.) C LP{d^), where p = 2{n - 
l)/(n — 2) if n > 3 or p < oo when n = 2. Thus the integral on the right-hand 
side will be well-defined if we have /iv in L2("-i)/"(5n) when n > 3 or LP{dn) 
for any p> 1 when n = 2. 

We say that a bounded, connected open set is a Lipschitz domain with 
Lipschitz constant M if the boundary is locally the graph of a Lipschitz func- 
tion. To make this precise, we define a coordinate cylinder Zr{x) to be a set of 
the form Zr{x) = {y : \y' — x'\ < r, |y„ — Xn\ < (1 -|- M)r}. We use coordinates 
{x', Xn) = (xi, x", Xn) G R X R"^^ X R and assume that this coordinate system 
is a translation and rotation of the standard coordinates. For each x in the 
boundary we assume that we may find a coordinate cylinder and a Lipschitz 
function (j) with Lipschitz constant M so that 

n n Zrix) = {{y', yn) ■■ vn > Hv')} n Zrix) 

dnnZrix) = :y„ = 0(y')}nZ^(a;). 

In the mixed problem, the boundary between D and N is another important 
feature of the domain. We assume that D is a relatively open subset of dil. and 
let A be the boundary (relative to d^l) of D. For each x in A, we require that 
a coordinate cylinder centered at x have some additional properties. We ask 
that there be a coordinate system {xi,x" ,Xn), a coordinate cylinder Zr{x), a 
function (j) as above and also a function ip : R""^ R with HVV'lloo < M so 
that 

Zr{x) nD = {{yi,y", y„) : yi > i^{y"), yn = Hv')} ^ Zr{x) 
Zr{x)nN = {{yi,y",yn) :yi<ip{y"), yn = (p{y')}nZr{x). 

We fix a covering of the boundary by coordinate cylinders {Zri{xi)}^i so that 
each ZiQQn (xi) is also a coordinate cylinder and let = minjrj : i = 1, . . . , N} 
be the smallest radius in the collection. We will use S{y) = dist(?/. A) to denote 
the distance from a point y to A. We will let Br{x) = {y :\y — x\ < r} denote 
the standard ball in R" and then Ar(x) = Br{x) D di} will denote a surface 
ball. Throughout this paper we will need to be careful of several points. The 
surface balls may not be connected and we will use the notation Ar{x) where 
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X may not be on the boundary. We use "^rix) to stand for Br{x) n $7. Since A 
is a Lipschitz graph, it has the property 

If X e A and < r < ro, then cr(Ar(x) T] D) > Af^V"-^ (2.2) 

Here and throughout this paper, we use a for surface measure. 

Our main tool for estimating solutions will be the non-tangential maximal 
function. For a > and x G (90, we define a non-tangential approach region 

by 

r(x) = {y edVL:\x-y\<{l + a) dist(y, dVL)]. 
Given a function u defined on $7, we define the non-tangential maximal function 

by 

v*{x) = sup X G dVL. 

It is well-known that for different values of a, the non-tangential maximal func- 
tions have comparable L^'-norms. Thus, the dependence on a is not important 
for our purposes and we suppress the value of a in our notation. 

Many of our estimates will be of a local, scale invariant form and hold for r 
less than a multiple of ro and with a constant that depends only on M and the 
dimension, n. Global estimates will also depend on the collection of coordinate 
cylinders which cover dil and the constant in (j2.ip . 

Before stating the main theorem, we recall the definition of atoms and 
atomic Hardy spaces. We say that a is an atom for the boundary if a is 
supported in a surface ball Ar{x) for some x in d^l, \\a\\L°°{dn) ^ ^/(^{^r{x)) 
and Jq^ ada = 0. 

When we consider the mixed problem, we will want to consider atoms for 
the subset A^. We say that a is an atom for A^ if a is the restriction to A^ of 
a function A which is an atom for (90. For a subset of A^ C 90, the Hardy 
space H^{N) is the collection of functions / which can be represented as ^ Xjaj 
where each aj is an atom for A^ and the coefficients satisfy X] l-^jl < This 
includes, of course, the case where A^ = 90 and then we obtain the standard 
definition. It is easy to see that the Hardy space H^(N) is the restriction to 
A^ of elements of the Hardy space H^{dQ). We give a similar definition for the 
Hardy-Sobolev space H^'^ of functions which have one derivative in H^. We 
say A is an atom for H^'^{dQ,) if A is supported in a surface ball Ar{x) and 
II Vtvl||/^oo(gQ) < l/a{Ar{x)). We say that A is an atom for H^'^ (D) ii A is the 
restriction to D of an atom for (90. Again, the space H^'^{D) is the collection 
generated by taking sums of atoms with coefficients in See the article of 
Coifman and Weiss [7] for more information about Hardy spaces. 

We are now ready to state our main theorem. 

Theorem 2.3 Let O, A^ and D be as described above. 

a) There exists qo > 1 so that for p in the range 1 < p < qo, given f^ S 
LP(A^) and fo G W^^P{D) we may find a solution u to (DUP with (Vu)* £ 
L^(90). The solution u satisfies 

\\iVu)*\\Lp(^QQ) < C{\\fN\\LP{N) + ll/D||iyi.P(D))- 
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The solution is unique in the class of functions with (Vn)* G LP{dVt). 

b) If fw lies in the Hardy space H^{N) and fo lies in the Hardy- Soholev 
space H^'^{D), then we may find a solution of the mixed problem u which 
satisfies 

||(Vn)*||2,i(an) < C{\\fN\\m{N) + \\fD\\m,-^{D))- 
The solution is unique in the class of functions with (Vu)* G L^{dQ). 

Proof. We begin by recalling that for the Dirichlet problem with data from 
a Sobolev space, we obtain non-tangential maximal function estimates for the 
gradient. This is treated for p = 2 hj Jerison and Kenig [15] and for 1 < p < 2 
by Verchota [331 El] • The Hardy space problem was studied by Dahlberg and 
Kenig [9]. Using these results, it suffices to prove Theorem l2.3l in the case when 
the Dirichlet data is zero. 

The existence for Neumann data from the atomic Hardy space and zero 
Dirichlet data is given in Theorem 14.141 and for data appears in section [71 
It suffices to establish uniqueness when p = 1 and this is treated in Theorem 



3 Higher integrability of the gradient of a 
weak solution. 

It is well-known that one can obtain higher integrability of the gradient of weak 
solutions. An early result of this type is due to Meyers [23]. We use the reverse 
Holder technique introduced by Gehring [T3] and, in particular, the formulation 
from Giaquinta [141 P- 122]. At a few points of the proof it will be simpler if 
we are working in a coordinate cylinder Z where we have that 90 n Z lies 
in a hyperplane. Thus, we will establish results for divergence form elliptic 
operators with bounded measurable coefficients as this class is preserved by a 
change of variable that will flatten part of the boundary. 

We define an operator P which takes functions on dfl to functions in by 

P/(x)=sup^/ I/I da 

and a local version of P by 

Prf{x)= sup / \f\da. 

r>s>0 JAs{x) 

On the boundary, we have that Pf is the Hardy-Littlewood maximal function 

Mf{x) = Pf{x) =snp^ [ \f\da, x G Oft. 
The following result is probably well-known, but we could not find a reference. 
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Lemma 3.1 For p > 1, 1 < q < pn/[n — 1), x G dVt and r < rg, we have 

iPrfl^dy] <c(^ [ fda^ . (3.2) 

The constant in this estimate depends only on the Lipschitz constant M and 
the dimension. 

Proof. We begin by considering the case where Q = {{y',yn) '■ Vn > 0} is a 
half-space. We use coordinates y = {y',yn) and we claim that 

Pf{y',yn) < Mf{y',0) (3.3) 

Pf{y) < c||/llLP(ac)yi'-"^/^ yn>o. (3.4) 

The estimate (I3.3P follows easily since As((y',y„)) C As((y',0)). To establish 
the second estimate, we observe that ii s < yn, then As{y) = and hence 



P/(y) = sup^/ \f\da<Cy 
s>y,i S JAsiy) 



(y) 



The second inequality follows from Holder's inequality. 

We claim that we have the following weak- type estimate for Pf, 

\{xGn: Pf{x) > A}| < C||/||^,(^^)A-W(-i), A > 0. (3.5) 

To prove (j3.5p . we may assume ||/||Lp(afi) = 1- With this normalization, the 
observation ([331) implies that {y' : Pf{y',yn) > A} = if y.„ > cA-^/^""^). 
Thus, we may use Fubini's theorem to write 

|{P/>A}| = / a{{y':Pf{y\yn)>X})dyn 
Jo 

< C / a{{y' :Mf{y',0)>cX})dyn 

Jo 

= CA^^"/^"""*"-* 

where we used ()3.3p . the weak-type {p,p) inequality for the maximal operator 
on R"~^ and our normalization of the L^-norm of /. 

From the weak-type estimate (|3.5p and the Marcinkiewicz interpolation the- 
orem we obtain 

< II/IIlp(R"-1)- (3-6) 

To obtain the estimate ()3.2p . we observe that if y G Br{x) then Br{y) C B2r{x) 
and hence 

Prf{y) < Pr{xA,A^)f){y), y e Br{x). 

Thus the result (13. 2p with 0, a half-space follows from ()3.6p and Holder's in- 
equality. 

Finally, to obtain the local result on a general Lipschitz domain, one may 
change variables so that the boundary is flat near Ar{x). I 
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We recall several versions of the Poincare and Sobolev inequalities. 

Lemma 3.7 Let il. be a convex domain of diameter d. Suppose that S <Z Cl is 
a set with the properties that for some r with < d < r we have a{S r]Br{x)) = 
r^~^ and a{S H Bt{x)) < Mt'^~^ for t > 0. Then for 1 < p < n, we have 



The constant depends only on M . 

Proof. We follow the proof of Corollary 8.2.2 in the book of Adams and Hed- 
berg pj, except that we have substituted the Riesz capacity for the Bessel ca- 
pacity in order to obtain a scale-invariant estimate. Following their arguments, 
we obtain that if u vanishes on S, then 

d^ 

Hx)\ < —(I,{\Vu\){x) + ||V^/|Up(f,)||/i(^)||^,,(^)). (3.8) 

Here Ii{f){x) = f{y)\x — dy is the first order fraction integral and fj, is 
any measure on 5 normalized so that /i(5') = 1. To estimate ||/i(/i)||^p'^^^ we 
use Theorem 4.5.4 of Adams and Hedberg [1] which gives that 

/ {h{ls)Y' dy<C j Wi:^^dfx 
where Wfp(x) is the Wolff potential of fi defined by 

roo 

{l^{B{x,t))tP-")^/^P-^^ dt/t. 



Our assumptions imply that with /i = r^^^o" denoting normalized surface mea- 
sure on S, we have /i(/u)(x) < C'r*^^"")/^^"^) where C depends only on M. 
Using this, the Lemma follows from ()3.8p . I 



The next inequality is also taken from Adams and Hedberg [H Corollary 
8.1.4]. Let 1/g + 1/n < 1 and assume that 17 is a convex domain of diameter 
d. We let u = j^udy and then we have 



(n+q)/n 



\u - u\'^ dy < C^y^ i^J I Vur«/('^+9) dyj . (3.9) 

Finally, we suppose that is a domain and ^rix) lies in a coordinate 
cylinder Z so that dO, H Z lies in a hyperplane and let u = u dy. Provided 

^r(a^) C Z, we have 

( [ \u-u\Ua] <c( [ \Vu\Pdy\ . (3.10) 

\JAr{x) J \J^r{x) J 

In this inequality, p and q are related by 1/q = 1/p — (1 — l/p)/{n — 1) and 
p>l. 
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Lemma 3.11 Let il., N, D be a Lipschitz domain for the mixed problem as 
defined above, suppose that It 2. ^|) holds, and let < r < rQ. Let u be a weak 
solution of the mixed problem for a divergence form elliptic operator with zero 
Dirichlet data and Neumann data f^. We have the estimate 



1/2 



\Vu\^dy 



< C 



{ 



I Vti| dy + 



1 



>.?i— 1 



i/p 



flida 



Here, p = 2 if n = 2 and p = 2{n — l)/(n — 2) for n > 3. The constant C 
depends only on M and the dimension n. 



Proof. Changing variables to flatten the boundary of a Lipschitz domain pre- 
serves the class of elliptic operators with bounded measurable coefficients, thus 
it suffices to consider the case where the ball Br{x) lies in a hyperplane. We 
may rescale to set r = 1. We claim that for 1/2 < s < t < 1, we have 



1/2 



\Vu\^dy 



< 



C 



it 



1/g / „ \ i/p 



\Vu\''dy 



+ 



ffj da 



(3.12) 

where we may choose the exponents p = 2{n — l)/(n — 2) and q = 2n/ (2n + 2) 
if n > 3 or p = 2 and q = 4/3 if n = 2. 

We give the details when n = 3. In the argument that follows, let e = 
(t — s)/2 and choose ?? to be a cut-off function which is one on Bs{x), supported 
in and satisfies \Vr]\ < C/e. We let v = rf{u — E) where E \s a, 

constant. If we choose E so that v G VF^'^(Q), the weak formulation of the 
mixed problem and Holder's inequality gives for 1 < p < oo 



Vu\'^rj^dy < C 



I \u-E\^\Vrf dy^\ [ \u-E\P'da] 
Jn \JAs+e(x) J 



+ 



As+,{x) 



\ 2/p' 

ffjda] 



(3.13) 



We consider two cases: a) Bs+e{x) r\ D = % and b) Bs+e{x) Ci D (/}. In 
case a) we may chooose E = u = ^ udy. We use the Poincare-Sobolev 
inequality (j3.9p and the inequality (j3.10p to estimate the first two terms on the 
right of (I3.13P and conclude that 



\Vu\^dy < C 



n + 2 



1 



\Vu\ "+2 dy 



+ 



np 

I Vu| dy 



2(np-n + l) 
pn 



+ 



flidcr 



Ai(x) 
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If n > 3, we may choose p = 2{n — 1) / (n — 2) and then we have that np/ {np — 
n + 1) = 2n/(n + 2) to obtain the claim. 

We now turn to case b). Since i?s+e(x) meets the set D, we cannot subtract 

1 2 

a constant from u and remain in the space of test functions, Wj^ (17). Thus, 
we let £' = in (I3.13p . We let u be the average value of u on ^ s+2e{x) and 
obtain 



v?\V7]\''dy<% 



\u — -up dy + 



Since i?s+e(x) n 7^ 0, our assumption on the boundary (j2.2p implies that we 
may find a point x so that B^{x) C Bt{x) and so that a{B^{x)r\D) > M~^e^~^. 
Using ()3.9p and the Poincare inequality of Lemma 13.71 we conclude that 



/ 



u^\Vri\^dy < C 



1 



(n+2)/n 



2-2ra/<j 



\Vu\ 



2n/{n+2) 



2/g 



\\/u\'^dy 



(3.14) 



for 1 < g < n. A similar argument using (j3.10p and Lemma 13.71 gives us 



i/p' 



\u\P da 



a + 2. 



Ax) 



< 



< c 



1/p' 



u — u 



P da 



+ u 



2 Ax) 



{np—n+l)/{np) 



+e 



'S's + 2e(x) 

1-n/q 



\^^\np/(np-n+l) ^ 
\ 1/? 



Ax) 



\Vu\'^dy 



(3.15) 



where the use of Lemma 13.71 requires that we have 1 < q < n. We use (j3.14p 
and (IXT51) in (f3T3]l and choose q = 2n/(n + 2) and p = 2{n - 2)/{n - 1) if 
n > 3. Once we recall that i — s = 2e, we obtain ()3.12p . 

Finally, we observe that the claim (I3.12P implies the estimate 



1/2 



\^u\ dy 



< C 



'^l/2(x) 



|Vn| dy + 



'fi(x) 



i/p 



fP 

Jn 



Ai(x) 



N 



(3.16) 



with p as in (j3.12p . This follows from an argument of Dahlberg and Kenig 
that can be found in Fabes and Stroock [121 pp. 1004-5]. The only change in 
their argument is to keep track of the extra term (J^^^^,) f^dyY^^ which we 

denote by A. To give the details, suppose that we have |Vn| dy = 1 and 

/ \l/q , , , , 

define I{t) = ( f^^^^^ | Vn|'? dy] where g < 2 is as in (fXT2|l . From (f3l^ . our 
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normalization in and Holder's inequality, we conclude that with = 2 — 2/q 
we have 

If I{t) < A for some t in [1/2, 1], then we have 1(1/2) < A. If I{t) > A for all t, 
then we have /(s) < CI{tf/{t - sf". Taking the logarithm of this expression, 
setting s = and integrating as in Fabes and Stroock [12j implies /(1/2) < C. 

When the dimension n = 2, the exponent 2n/(n + 2) is 1 and it is not clear 
that we have (|3.9p as used to obtain (I3.14p . However, from (|3.9p and Holder's 
inequality we can show 



1/2 



\u — u\ dy 



< C 



\ 3/4 



This may be substituted for (|3.9p in the above argument to obtain (j3.12p when 
n = 2. I 

Lemma 3.17 Let u he a weak solution of the mixed problem with zero 

Dirichlet data and Neumann data f in U'{N) which is supported in N riAr{x) 
with r < rQ. There exists po > 2 so that for t in [2,po); we have the estimate 



i/t 



ix) 



\Vufdy 



< C 



j 



|Vn| dy + 



n/{t{n-l))- 



ft{n-l)/n 



Acr{x) 



The constant in this estimate depends on M and n. 

Proof. According to Lemma [3. IH Vu satisfies a reverse Holder inequality and 
thus we may apply a result of Giaquinta [14', p. 122] to conclude that there 
exists po > 2 so that we have 



4 jVuI'dy 



1/t 



< C 



j 



\Vu\ dx + 



/ {pfpY/f 



l/p 



dy 



for t in [2,po) and p as in lemma [3.111 From this, we may use Lemma l3.ll to 
conclude the estimate of the Lemma. I 



4 Estimates for solutions with atomic data 

We establish an estimate for the solution of the mixed problem when the Neu- 
mann data is an atom and the Dirichlet data is zero. The key step is to establish 
decay of the solution as we move away from the support of an atom. We will 
measure the decay by taking L'^-norms in dyadic rings around the support of 
the atom. Thus, given a surface ball Ar(x), we define = iS.2k.i,{x)\/S.2k-i^{x) 
and define = ^2'=r(^) \ ^2*-ir(^) = ^^t,k-i^3- 
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Theorem 4.1 Let u be a weak solution of the mixed problem with Neu- 

mann data fN = aan atom for N which is supported in Ar{x) and zero Dirich- 
let data. For 1 < q < pq/2, we have the following estimates 

\Slu\ida\ <Ca[^Ar{x))-^l'<! (4.2) 



and for k > 3 



\ 1/9 



Q' \Vu\'ida^ % C2-"V(Sfe)-W. (4.3) 



Here, a is as Lemma \4. 7[ po is the exponent from Lemma 3.17, and the constant 
C depends on 

If r < ro and x is in dQ, then we may construct a star-shaped Lipschitz 
domain r2j.(x) = Z^.(x)nil where Zr{x) is the coordinate cyhnder defined above. 
Given a function u defined in fi, x G dil., and r > 0, we define a truncated non- 
tangential maximal function v* by 

v*(x) = sup \v{y)\. 

y&{x)nBr{x) 

Lemma 4.4 Suppose that x € dfl and r < tq. Let u be a harmonic function 
in Q,4r{x). If Vu G L^(ri4r(x)) and du/dv is in L'^{dQ. R dQ^rix)), then we 
have Vu E L'^{Ar{x)) and 



\Vu\^dy 



[ {{Vu);fda<cif 
JAr(x) \JdS 





du 




I 




da + - [ 


\J anndfiirix) 


r Jn 



Proof. Since the estimate only involves Vn, we may subtract a constant from 
u so that j^^j udy = 0. We pick a smooth cut-off function r] which is one 

on Z-ir{x) and zero outside Z^j-ix). Since we assume that Vn is in L'^{Q), it 
follows that A(r]u) = 7]Au + 2Vu ■ Vrj is in L^(ri). Thus, with H the usual 
fundamental solution of the Laplacian, w = 3* {A{r]u)) will be in the Sobolev 
space VK^'-^(R"). Next, we let v be the solution of the Neumann problem 

in 174,. (x) 
on 504r(x). 

According to Jerison and Kenig [15], the solution v will have non-tangential 
maximal function in L'^{dQ). By uniqueness of weak solutions to the Neumann 
problem, we may add a constant to v so that we have rju = v + w. As w and all 
its derivatives are bounded in Q2rix) and the non-tangential maximal function 
of Vv is in L^(9il4^(x)), we obtain the Lemma. I 
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In the next lemma and below, we let V^n = Vn — z^Vn • v denote the 
tangential component of the gradient. The proof of this lemma is identical to 
the proof of Lemma [ 



Lemma 4.5 Suppose that x G 90 and r < tq. Let u be a harmonic function 
in ^2r{3:)- If^u G L'^{^4r{x)) and Vtu is in L?'{dQ. Pi dil^rix)), then we have 
Vu e L2(Ar(x)) and 

I {{Vu);f da < C I [ \Vtu\'^da + - [ |Vupdy). 

JArix) \JdQr\dn4r{x) ^ Jn4r{x) J 

The following weighted estimate will be an intermediate step towards our 
estimates for solutions with atomic data. In the next lemma is a bounded 
Lipschitz domain and the boundary is written di} = D L) N. Recall that 6{x) 
denotes the distance from x to the set A. 

Lemma 4.6 Let e G R, r < tq, and suppose that u is a weak solution of the 
mixed problem with /z? = and /at in L'^(N). There is a constant C so that 
the solution u satisfies 

[ {{^u)lsf6^''da<C\[ \fNfS^''d(j+ f dy] . 

for constants c and C which depend only on n and M . 

The next lemma uses a Whitney decomposition and thus it is simpler if we 
use surface cubes, rather than the surface balls used elsewhere. A surface cube 
is the image of a cube in R""^ under the mapping x' (x', (t){x')). Obviously, 
each cube will lie in a coordinate cylinder. 

Proof. Note that if 5{y) > r for all y in Ar(x), then this Lemma follows 
directly from Lemma 14.41 or Lemma 14.51 The Lemma is more interesting when 
Ar(x) is close to A. To prove the Lemma, we make a Whitney decomposition 
of 90 \ A. Thus we write dVL = A U (UjQj). The collection of surface cubes 
{Qj} has the following properties: 1) for each j, we have either Qj C D or 
Qj C (A^ \ A), 2) we have c6{x) < diam((5) < (5(x)/4 for each x in Q, 3) if we 
define T{Q) = {x £ Cl : dist{x,Q) < diam((5)}, then the sets {T{Qj)} have 
bounded overlaps and thus 

J2xTiQ,)<C{n,M). 
We let rj = diam(Qj) and use the estimates of Lemma [4. 51 and [4. 41 to obtain 

/ iVu)*.\Vu\'^da <C [ 1^ da + - [ \Vu\^ dy 
jQj ' \j2Q,nN ou rj Jt{Qj) 

Now we multiply by r^*^, use that rj ~ 5{x) for x £ T{Qj), sum on j and use 
that the family {T{Qj)} has bounded overlaps to obtain the Lemma. I 
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An important part of the proof of our estimate for the mixed problem is 
to show that a solution with Neumann data an atom will decay as we move 
away from the support of the atom. This decay is encoded in estimates for the 
Green function for the mixed problem. These estimates rely in large part on 
the work of de Giorgi [lOj, Moser [25] and Nash [26] on Holder continuity of 
weak solutions of elliptic equations with bounded and measurable coefficients, 
and the work of Liftman, Stampacchia and Weinberger [T9j who constructed 
the fundamental solution of such operators. Also see Kenig and Ni [17J for the 
construction of a fundamental solution in two dimensions. Given the free space 
fundamental solution, the Green function may be constructed by reflection in a 
manner similar to the construction given for graph domains in [18]. A similar 
argument was used by Dahlberg and Kenig ^ in their study of the Neumann 
problem. 

Lemma 4.7 We consider the mixed problem in a Lipschitz domain with D and 
N satisfying our standard hypotheses. There exists a Green function G{x, y) for 
which satisfies: 1) If Gx{y) = G{x,y), then Gx is in w]f' {^\ Bj.{x)) for 

all r > 0, 2) AGx = 6x, the Dirac 6 -measure at x, 3) If fN He-s in Wj^ ' (dO,), 
then the solution of the mixed problem with fD = can be represented by 

u{x) = -{fN,Gx)dn, 

4) The Green function is Holder continuous away from the pole and satisfies 
the estimates 

\Gix, y) - Gix, y')\ < ^^^^1", , \x - y\ > 2\y - y'\, 

\Gix,y)\<- n>3, 

\x — y\"' ^ 

and 

|G(x,y)| <C(l + |log|x-y||), n = 2. 

Lemma 4.8 Let u be a weak solution of the mixed problem with Neumann data 
f in L^{N) where p = (2n — 2)/n for n > 3. Then we have the estimate 

\Vu\Uy<C\\f\\l,^j,y 
\yu\'dy<C\\f\\l,^^y 



If n = 2, we have 



In each case, the constant C depends on $7. 

Proof When n > 3, we use that W^^'^''^{dn) C L2("-i)/("-2) (^j^)^ Transpos- 
ing to the dual tells us that L2("~^)/"(5il) C W^^^'^'^{dQ,) and since the weak 
solution of the mixed problem satisfies 



n 



\Vu\^dy<C\\f\\] 



-1/2,2 



(an) 
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the Lemma follows. 

When n = 2, the proof is the same except that we do not have that 
C L°°{dQ). However, we do have the embedding W^^'^''^{dQ) C 
BMO{dO,) and this gives the conclusion of the Lemma for n = 2. I 



Proof of Theorem \4-i\ We begin by considering 



\Vu\Ua 



to obtain an estimate for the gradient of the solution near the support of the 
atom. To begin we let g > 1 and use Holder's inequality to estimate 



\ 1/9 



1/2 



\Vu\^5^-'da 



< Cr 



6 2-q da 



2-q 
2q 



1/2 



Asr{x) 



This requires that q and e satisfy q{e — l)/(2 — q) > —1 or q < 1/(1 — e/2). 
Next we use Lemma 14.61 to bound the weighted {6^^'' da)-norui of Vu. This 
gives us 



1/9 



\Vu\ida 



< C 



Asr{x) 



da 



Ar{x) 



1/2 



+ 



<i>Crix) 



\ 1/2 

\Vu\'^5-'dy] \r 



("-l)(i-5)+^ 



To estimate the solid integral in this last expression, we use Holder's inequality 
again and find 



1/2 



Vu\'^5-'dy\ < 



•^Cr{x) 



1/p 



^ '"^Crix) 



\Vu\Pdy 



^-ep/(p-2) 



^Cr{x) 



\Vu\Pdy 



i/p 



^Ct{x) 



The last inequality holds when ep/{p — 2) <2orp>2/(l — e/2). Now we may 
use the previous two displayed inequalities, Holder's inequality and Lemma [3.17l 
to obtain 



1 



\ 1/9 

|Vn|'?d(T <C 



1 



, ,2 

|Vn| dy \ + r 



1— n 
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We have used the normaUzation of the atom to estimate the term involving 
the Neumann data in Lemma 13.171 Finally, we may use Lemma 14.81 and the 
normalization of the atom o to obtain the estimate (r"" J | Vup dy^/"^ < Cr^~^ 
which gives the estimate (14. 2p . Examining the conditions onp > 2/(1 — e/2) and 
(7 < 1/(1 — e/2) we see that we obtain this estimate precisely when 1 < q < pq/2 
and pq is the exponent from the reverse Holder argument in Lemma 13.171 

Our next goal is to estimate /^^^ \Vu\'^ da for k sufficiently large. If 1 < g < 2 
and g < 1/(1 — e/2), we may use Holder's inequality to obtain 

Vul^daJ < \Vu\^5^-'d(jj yj d^^dcjj 

< c(^j |Vn|25i-^day^'(2V)("-i)(i-5)+(-i)/(i9) 

This requires that (e — l)q/{2 — q) > —1. Next, Lemma 14.61 gives that 

J \Vu\'^6^-'daj <Ci^J \Vu\^6-'dyj . (4.10) 

To estimate the integral over fJ, we use Holder's inequality and then the reverse 
Holder estimate in Lemma 13.171 to obtain 

\Vu\'^5-'dy\ < C( \Vu\Pdy] (2'^r)"(p"^)"^ 

1/2 



< c(^j \Vu\'^dy^ (2'=r)-f. (4.11) 



'Sk 

Here we have used that {f^^ 5-^p/(p-2) dyy/P < C(2'=r)"(i/P-i/2)-€/2 provided 
p > 2/(1 — e/2). We appeal to the Caccioppoli inequality, recall that the mixed 
data for u is zero outside of Ar{x), and obtain that 

N 1/2 . . ^ X 1/2 

^^\Vu\^dyj <C^(^Jju\^dyj . (4.12) 
From Lemma 14.7^ we represent u by 



uiy) = - G{y,z)a{z)dz. 

JNnArix) 

If Ar{x) intersects D, we have that G is zero at some point in Ar{x). Thus, the 
normalization of a, ||a|| ^1(9^2) = 1, and the Holder continuity of G in Lemma 
14.71 part 4) imply that 

|n(y)| < C7r"/|x-yr-2+", if \x - y\ > 2r. (4.13) 

If Ar{x) C N, then we may use that a has mean value zero to write 



u{y)=\ {G{y,z) - G{y,x))a{z)dz 

JNnArix) 
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and obtain the estimate (j4.13p from part 4) of Lemma 14. 7i Using the estimate 
(j4.13p in (|4.12jl and the estimates (l4.9H4.lll) gives the result of the theorem. 
We note that if po is the exponent for the reverse Holder inequality then the 
conditions on e allow us to again obtain results for 1 < q < pq/2. I 

We now show that the non-tangential maximal function of our weak solu- 
tions lies in when the Neumann data is an atom. 

Theorem 4.14 Let fj^ he in H^{N), then there exists u a solution of the 
mixed problem with Neumann data /at and zero Dirichlet data and this solution 
satisfies 

ll(Vn)*||ii(af,)<C||/||Hi(^). 

Proof. We begin by considering the case when /jv is an atom and we let u be 
the weak solution of the mixed problem with Neumann data an atom a and 
zero Dirichlet data. The result for data in (N) follows easily from the result 
for an atom. 

We establish a representation for the gradient of u in terms of the boundary 
values of u. Let x E O and j be an index ranging from 1 to n. We claim 

+^i^--)Pda. (4.15) 
oyj ov 

If u is smooth up to the boundary, the proof is a straightforward application 
of the divergence theorem. However, it takes a bit more work to establish this 
result when we only have that u is a weak solution. 

Thus, we suppose that r/ is a smooth function which is zero in a neighbor- 
hood of A and supported in a coordinate cylinder. Using the coordinate system 
for our coordinate cylinder, we set Uriy) = u{y + Ten)- Applying the divergence 
theorem gives 

- / r^{^{x - ■)-^ - VH(x - •) • VurVj + (x - .)^) da 
Jdn OL' dyj dyj du 

du f du 

= Vix)^{x)- / Vt? • VH(x - •)— ^ 
CXj Jq oyj 

-VyE{x - •) • Vurp- 
dyj 

dE 

+—(x--)Vur-Vr]dy. (4.16) 
9yj 

Thanks to the truncated maximal function estimate in Lemma l4.6t we may let 
r tend to zero from above and conclude that the same identity holds with Ur 
replaced by u. Next, we suppose that rj is of the form r/(/>e where (/>e = on 
{x : 6{x) < e} and (pe{x) = 1 on {x : 5{x) > 2e} and we have the estimates 
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|V(/>e(2;)| <C/e. Since we assume the boundary between D and is a Lipschitz 
surface, we have the following estimate for e sufficiently small 

\{x : 6{x) < 2e}| < C7e^ (4.17) 

Using our estimate for V(/>e and the estimate (I4.17p . we have 

\ [ vVcj),-VE{x--)-^dy\<C I [ \Vu\^dy] 

Jn Oyj \J {y:5{y)<2e} J 

and the last term tends to zero with e since the gradient of a weak solution lies 
in L^(0). Using this and similar estimates for the other terms in (|4.16p . gives 

hm - [ Vi^^ii) ■ VS(x - - VyE{x - •) • Vn^ 

e^o+ 7n oyj dyj 

d'^ f du 

+—{x--)Vu-V{(t>ev)dy = - V{(P,r])■VE{x-■) — 
^yj Jn dyj 

yj-f \ V7 ^^^^ 
— \/=l(x — •) • Vn— — 

dE 

^dy~^^ - ■)'^u ■ Vrjjdy. 

Thus we obtain the identity ()4.16p with Ur replaced by u and without the 
support restriction on ij. Finally, we choose a partition of unity which consists 
of functions that are either supported in a coordinate cylinder, or whose support 
does not intersect the boundary of fi. Summing as r] runs over this partition 
gives us the representation formula (j4.15p for u. As we have Vit G L'^{dVt), 
it follows from the theorem of Coifman, Mcintosh and Meyer [6] that (V«)* 
lies in L'^{dQ,). However, a bit more work is needed to obtain the correct 
estimate for (Vu)*. 



We claim 



Ian 
du du 



du , 

da = 



/ J^jT, Vi^da = 0. 

Jan oVi oyj 



Since (Vu)* lies in L'^{dQ), the proof of these two identities is a standard 
application of the divergence theorem. Using these results and the estimates 
for Vu in Theorem 14. 11 allows us to show that du/du and Ujdu/dyi — uidu/dyj 
are molecules on the boundary and hence it follows from the representation 
formula (j4.15p that (Vu)* lies in L^{d^) and satisfies the estimate 

\\{^nr\\LU9n)<C. 
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5 Uniqueness of solutions 



In this section we establish uniqueness of solutions to the mixed problem (jl.lh . 
We use the existence result established in section d] and argue by duality that 
if n is a solution of the mixed problem with zero data, then u is also a solution 
of the Dirichlet problem with zero data and hence is zero. 

Theorem 5.1 Suppose that u solves the mixed problem ( li. j|) with data /at = 
and fn = 0. If (Vu)* G L^{dn), then n = 0. 

We recall that G. Verchota constructed a sequence of approximating do- 
mains in his dissertation Theorem 1.12] j33l Appendix A]). We will need 
this approximation scheme and a few extensions. Given a Lipschitz domain 
Q, Verchota constructs a family of smooth domains {0,^} with 0^ C $7. In 
addition, he finds bi-Lipschitz homeomorphisms A^ : di^k — *■ dQ which are 
constructed as follows. 

We choose a smooth vector field V so that for some 5 = 5{M), V ■ u < —5 
a.e. on dVl and define a flow for t small by ^f{x,t) = V{f{x,t)), f{x,0) = x. 
One may find ^ > so that 

O = {f{x, t) : X e dn,-^ <t < ^} (5.2) 

is an open set and the map {x,t) — > f{x,t) from dO, x — > O is bi- 
Lipschitz. Since the vector field V is smooth, we have 

Df{x,t) = ln + 0{t) (5.3) 

where In is the n x n identity matrix and DF denotes the derivative of a 
map F. In addition, we have a Lipschitz function tk{x) defined on di} so that 
^kix) = f{x, tk{x)) is a bi-Lipschitz homeomorphism A^ : dQ dH.^- We have 
a collection of coordinate cylinders {Zi} so that each Zi serves as a coordinate 
cylinder for dil. and for each of the approximating domains d^l^- If we fix a 
coordinate cylinder Z, we have functions <j) and (p^ so that dQCiZ = {{x' , 4>{x')) : 
x' G R"'i}nZ and n Z = {(x', 0fc(x')) : x' G W-^]r\Z. The functions (/>fc 
are C°° and || V'0A,.||/^oo(j^n-i) is bounded in k, limfe^oo V0fc(2;') = V'</>(x') a.e. 
and (pk converges to (j) uniformly. Here we are using V' to denote the gradient 
on R"~i. 

We let TT : R" — > R""^ be the projection 7r(x',x„) = x' and define Sk{x') = 
TT(Ak{x' , (p{x'))). According to Verchota, the map Sk is bi-Lipschitz and has a 
Jacobian which is bounded away from and cxd. We let Tk denote S^^ and 
assume that both are defined in a neighborhood of it{Z). We claim that 

lim DTk{Sk{x)) = In-i, a.e. in 7r(Z), (5.4) 

fc— >oo 

and the sequence ||Z)Tfe||ioo(^(^)) is bounded in k. 

To establish (j5.4p . it suffices to show that -05^ converges to In-i and 
that the Jacobian determinant of DSk is bounded away from zero and in- 
finity. The bound on the Jacobian is part of Verchota's construction (see 133\ 
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p. 119]). As a first step, we compute the derivatives of tk{x' , (p{x')). Since 
f{{x' , (l){x')),tk{x' ,(f){x'))) lies in OOfc, the derivative is tangent to dQk and we 
have 

{—f{x, (f){x'),tk{x', (pix')))) ■ Uk{y) = 0, a.e. in 7r(Z), 

UXi 

where y = {Sk{x'),4>k{Sk{x'))) and Vj. is the normal to 80,^- Solving this 

dxi ' 



equation for i^rtk gives 



^tk{x'A{x')) = -{V{y).Uk{y))-'{^{{x\<P{x'))Mx'A{x'))) 

OXi \OXi 

+ ^ix')-^{{x' ,4>ix')),tk{x' ,(j){x')))^ ■ Ukiy). 

Since liuik^ootk = uniformly in tt{Z), (j5.3p . and i^k{y) converges pointwise 
a.e. and boundedly to we obtain that 

d 

lim — — tk(x' ,4>(x')) = 0, a.e. in 7r(Z). (5.5) 

fc^oo OXi 

Computing 

-^Jiix',<Pix'))Mx',cf>ix'))) 

= ^(x', Hx'),tkix', Hx'))) + 4>{x')M^', 'A(x'))) 

+V{f{{x', 0(x')), tfc(x', </,fc(x'))))^tfc(x', <A(x')). 

OXi 

Given ([O]), ([53]), and recalling that Sk{x') = TT{f{x' , (j){x')),tk{x' , (p{x'))), ([531) 
follows. 

Lemma 5.6 // w is the weak solution of the mixed problem with Neumann 
data an atom and zero Dirichlet data, then we have 



an 



U—da < Cw\\u\\yyl,l(^gQ^y 



k 



Proof. This may be proven using generalized Riesz transforms as in |34t Sec- 
tion 5]. Verchota's argument uses square function estimates to show that the 
generalized Riesz transforms are bounded on L^. In the proof of this Lemma, 
we need that the Riesz transforms of w are bounded. From the estimate for 
the Green function in Lemma 14.71 and the representation ol w = —{G,a)Qfi, 
we conclude that w is Holder continuous. The Holder continuity, and hence 
boundedness, of the Riesz transforms of w follow from the following character- 
ization of Holder continuous harmonic functions. A harmonic function n in a 
Lipschitz domain is Holder continuous of exponent a, < a < 1, if and only 
if sup2,gQ dist(x, (90)^~°| Vn(x)| is finite. I 
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We will need the following technical lemma on approximation. The proof 
relies on the approximation scheme of Verchota outlined above. In our applica- 
tion, we are interested in studying functions in Sobolev spaces on the family of 
approximating domains. Working with derivatives makes the argument fairly 
intricate. 

Lemma 5.7 If u satisfies (Vu)* G and \/u has non-tangential limits 

a.e. on dO,, then we may find a sequence of Lips chitz functions Uj so that 

lim ll-u - Uj\\whi{ank) ^ ^'/j- 

If u\gn is supported in N , then C/j|an is also supported in N . 
The constant C may depend on Q.. 

Proof. To prove the Lemma, it suffices to consider a function u which is zero 
outside one of the coordinate cylinders Z as given in Verchota's approximation 
scheme. We have u{x' ,4>{x')) S VF^'-'^(R"~-'^). Hence, there exists a sequence of 
Lipschitz functions Uj so that /j^n-i | V'n(x', — Vuj{x' .,(j){x'))\ dx' < 1/j 

where V' denotes the gradient in R"^^. We extend uj to a neighborhood of 
dn by 

Uj{f{x,t)) = 7]{f{x,t))u,{x), xedn 

where r/ is a smooth cutoff function which is one d^l and supported in the set 
O defined in (lOI) . 
We consider 

/ \Vu{x', Mx')) - VU,{x', Mx'))\ dx' 

J-k{Z) 

< [ \V'u{x',^k{x')) -Vu{x',^{x'))\dx' 

+ / \Vu{x', - V'ujix', (/>(x'))| dx' 

Jtt(Z) 

+ / \Vuj{x' ,(l){x')) - VUj{x' ,^k{x'))\dx' 
= Ak + B + Ck. 

We have that Wvuk^^Ak = since we assume that (Vu)* G L^{d^), Vu has 
non-tangential limits a.e., and V'(/>fc converges pointwise and boundedly to V'(/). 
By our choice of Uj, we have B < C/j. Finally, our construction of Uj and our 
definition of imply that Uj{x' , (pkix')) = Uj{Tk{x'), (p{Tk{x'))) and hence we 
have 

Ck < [ \{In-i-DTk{x'))Vuj{x',^{x'))\dx' 

Jn{Z) 

+ / \Dn{x'){Vu.j{x', (t^{x'))) - v'u,{n{x'), HTkix'm dx' 

Jn{Z) 
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We have that hm, 



since Vu^ is bounded and (15.41) holds. Since 



Tk{x') converges uniformly to x', DTk is bounded, and the Jacobian of Sk is 
bounded, we have that linifc^oo C'fc,2 = 0. I 

Proof of Theorem \5.1[ We let u be a solution of (jl.ip with zero data /at and 
fn and we wish to show that u is zero. We fix a an atom for and let w he a 
solution of the mixed problem with Neumann data a and zero Dirichlet data. 
Our goal is to show that 

/ auda = 0. (5.8) 
Jn 

This implies that u is zero on and then Dahlberg and Kenig's result for 
uniqueness in the regularity problem [9j imply that n = in 

We turn to the proof of (|5.8p . Applying Green's second identity in one of 
the approximating domains Ofc gives us 



an, 



du , 
w— da 
ov 



dw , 

Oh' 



k = 1,2.... 



(5.9) 



We have (Vu)* is in L^{dO,) and w Holder continuous and hence bounded. 
Recalling that w is zero on D and du/dv is zero on N , we may use the dominated 
convergence theorem to obtain 



lim 



an, 



w— da = (J. 
ov 



(5.10) 



Thus, our claim will follow if we can show that 

dw 



lim 



dVi, 



u—— da 
ov 



ua da. 



an 



(5.11) 



Note that the existence of the limit in (|5.11|) follows from (j5.9|) and (|5.1U|) . We 
let Uj be the sequence of functions from Lemma 15.71 and consider 



ua da 



an 



lim 

fc— >oo 



f dw , 








< 


j 






Jan 



ua da 



lim 

fe— +cx> 



an. 



ov 



lim sup 

fc— >oo 



an^ 



ov 



Because we have that (Vw)* is in L^{dU) and Uj is bounded, we may take the 
limit of the first term on the right of (j5.12p and obtain 



ua da 



an 



lim 



an, 



dw 
u—— da 
ov 



< 



an 



{u — Uj)a da 



+ C/j<C/j. 



According to Lemmata 15.61 and 15.71 the second term on the right of (j5.12p is 
bounded by C^/j. As j is arbitrary, we obtain (j5.1ip and hence the Theorem. 
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6 A Reverse Holder inequality at the bound- 
ary 

In this section we establish an estimate in U'{d^) for the gradient of a solution 
to the mixed problem. This is the key estimate that is used in section [7] to 
establish L^-estimates for the mixed problem. 

Lemma 6.1 Let u he a solution of the mixed problem M-l\) with zero Dirichlet 
data and Neumann data Jn- Let po be as in Lemma \3.17\ and assume /at G 
^Po{n-i)/n(jV). Then for 1< q < po/2, x e ft and r < tq, we have 



1/9 



(/ {Vu);Ua] <C -f \Vu\dy+[-f 

\JAr{x) J J^Crix) \JNnAcr{x) 

The constant in this local estimate depends on n and M . 



fPo(n-l)/n 
J N 



{n-l)/(npo)' 



da 



Proof. By rescaling, we may assume r = 1. Fix x G 17 let A,. = Aj.(2;), 
Br = Br{x) and = Br n Q to simplify the notation. We let e satisfy 
e < 2 — 4/po where po is as in Lemma 13.171 Using Holder's inequality, we have 



)g/{2-g) 



l/g-l/2 



where the last integral will be finite if q < 1/(1 — e/2). From Lemma |4.6| we 
obtain 

r \l/2 / r N 1/2 / . n 1/2 

/ \Vu\'^S^-'da) <C{ f^5^-'da) + / S-'\Vu\'^dy] . 

JA2 / \JA4, ) \J^4 ) 

By Holder's inequality and Lemmata 13.171 and 13.111 we may find an exponent 
p > 2 so that 



lj>4 



1/2 



< 



< c 



\^Ip f f \ l/'i-n/p 

Vu\'Pdy] / 6-'P/^P~^Uy 

, , ^ in-^)/inp)' 

Vu\dy+[ I fl'^-'^'-da 



where our choice of e implies that the integral of 6 ^) is finite. Combining 

the previous two displayed equations, we obtain 



\ 1/9 

iViil^dcT <C 



\Vu\ dy 



f 



po{n-l)/n 
N 



da 



"/{("-l)Po)' 



(6.2) 

We now obtain an estimate for (Vn)* in Towards this end, let r] be 

a cut-off function which is one on the ball Bi and is supported in B2- As the 
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Neumann data for our solution u is an atom, we have that (Vu)* G L^[d^) 
from Theorem 14.11 Hence, we have the representation formula 

an^ydyi ^ '^'dyj dyi ^ dyi^^ yj ^ dyj^'J ^ 

= -(„u)(z) + / ^^^f, _ .) _ _.)^^ + ^(,_.)^^dy 

Jn dyi dyj dyi dyi dyi dyj yj dyj dyi 

This representation formula and the theorem of Coifman, Mcintosh and Meyer 
on the Cauchy integral [6j imply that 

ll(V^x):iU,(A,/,) < C(||V^x|U,(A,) + / |Vn| dy). 

The Lemma follows from this estimate and (16.21). I 



7 Estimates for solutions with data from L^, 

p> 1 

In this section, we use the following variant of an argument developed by Shen 
|29j to establish L^-estimates for elliptic problems in Lipschitz domains. This 
argument appeared earlier in work of Peral and Caffarelli [5j. 

As the argument depends on a Calderon-Zygmund decomposition into dyadic 
cubes, it will be stated using surface cubes rather than the surface balls Ar{x) 
used elsewhere in this paper. 

Let Qo be a cube in the boundary and let F be defined on AQq. Let the 
exponents p and q satisfy 1 < p < q. Assume that for each Q C Qo, we may 
find two functions Fq and Rq defined in 2Q such that for 



/ 

/ \Rq\ 

J2Q 



Fglda < C{4 \f\Pda 
2Q VJ4Q 

1/9 

"da] < C 



i/p 



F\ < \Fq\ + \RqI 

/ i/r 

/ \F\da+(J \f\PdaY' 



(7.1) 
(7.2) 

(7.3) 



Under these assumptions, for r in the interval {p,q), we have 



/ 

JQo 



1/r 

iF^da] <C 



/ \F\da+U i/r 



da 



l/r 



The constant in this estimate will depend on the Lipschitz constant of the 
domain and the constants in the estimates in the conditions (|7.2H7.3|) . The 
argument to obtain this conclusion is more or less the same as in Shen |291 
Theorem 3.2]. We are not able to use Shen's result directly as we have results 
in Hardy spaces rather than L^-spaces. 
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In our application, we will let 4Qo be a cube with sidelength comparable 
to tq. We let u he a solution of the mixed problem with Neumann data / in 
U'(N) and Dirichlet data zero. Since IJ'{N) is contained in the Hardy space 
H^{N), we may use Theorem 14.11 to obtain a solution of the mixed problem 
with Neumann data / and zero Dirichlet data. Let F = (Vn)* and given a cube 
Q C Qo and with diameter r, define Fq and Rq as follows. We let f^Q = 
if 4Q n D / and Uq = f^Q f da if 4Q C N. Set g = xaqU - Iaq) and 
h = f — g. As both g and h are elements of the Hardy space H^{N), we may 
solve the mixed problem with Neumann data g or h. We let v be the solution 
with Neumann data g and w be the solution with Neumann data h. According 
to the uniqueness result Theorem 15.11 we have u = v + w. We let Rq = (Vw)* 
and Fq = (S/v)* so that (|7.ip holds. We turn our attention to establishing 
(lEl and 

To establish ()7.2p . observe that the H^-norm of g satisfies the bound 

\\9\\hHn) < C\\f\\L.i4Q)cr{Q)'/P'. 

With this, the estimate (j7.2p follows from Theorem 14.11 Now we turn to the 
estimate (I7.3P for Fq = (Vw)*. We note that the Neumann data h is constant 
on 4(5 n A^. We define a maximal operator by taking the supremum over that 
part of the cone that is far from the boundary, 

(y'w)^{x) = sup \V'w{y)\ 

y&r{x)nBcr{x) 

where C is to be chosen. 

We first observe that we have 

{Vw)l{x) <C-f (Vw)* da, X G 2Q. (7.4) 

The estimate for (Vw)'^^ uses the local estimate for the mixed problem in 
Lemma l6. II to conclude that 




We choose the constant C in the definition of (Vw)^ large in order that (|7.5p 
holds. Recall that T{Q) was defined at the beginning of the proof of Lemma 
14.61 From the estimates (|7.4p and (j7.5p . we conclude that 
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We have (Vw)* < (V v)* + (V u)* and hence we may estimate the term involving 
the (Vw)* by 



/ {Vwyda<-f {Vuyda+-f (Vv)* da < -f {Vu)*da+c(i fP da] 

J\Q J\Q JA.Q JaQ xJaQ J 



Combining this with (j7.6|) gives (j7.3|) . 

Applying the technique of Shen outlined above gives the L^-estimate and 
thus we obtain the following theorem. 

Theorem 7.7 Suppose that 1 < p < pq/2 where pQ is as in Lemma 1 7\ Let 
Jn £ LP[N) and let u he the solution of the mixed problem with Neumann data 
/at and zero Dirichlet data. 
The solution satisfies 

ll(V'")*llLP(9n) < C\\fN\\LP{dn)- 

The constant depends on the number of coordinate cylinders and p as well as 
M. 



8 Further questions 



This work adds to our understanding of the mixed problem in Lipschitz do- 
mains. However, there are several avenues which are not yet explored. 

1. Can we study the inhomogeneous mixed problem and obtain results simi- 
lar to those of Fabes, Mendez and Mitrea [Tl] and M. Mitrea and I. Mitrea 



2. Is there an extension to p < 1 as the work of Brown [2\? 

3. Can we study the mixed problem for more general decompositions of the 
boundary, = DU N? To what extent is the condition that the bound- 
ary between D and be a Lipschitz graph needed? 

4. Can we extend these techniques to elliptic systems and higher order elliptic 
equations? 
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